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Mutually unbiased bases and discrete Wigner functions are closely, but not uniquely related. Such 
a connection becomes more interesting when the Hilbert space has a dimension that is a power of a 
prime N = d n , which describes a composite system of n qudits. Hence, entanglement naturally enters 
the picture. Although our results are general, we concentrate on the simplest nontrivial example of 
dimension N — 8 = 2 3 . It is shown that the number of fundamentally different Wigner functions 
is severely limited if one simultaneously imposes translational covariance and that the generating 
operators consist of rotations around two orthogonal axes, acting on the individual qubits only. 



PACS numbers: 



I. INTRODUCTION 



Quasiprobability distributions have been useful in 
quantum mechanics both as visualization tools and for 
computational purposes. Typical examples include the 
Glauber-Sudarshan HQ, the Husimi 0, and the Wi ancr 
function Q . These are all functions of continuous param- 
eters that map quantum states onto a continuous phase 
space. 

The growing interest in quantum information has fu- 
eled the search for discrete phase-space counterparts of 
these distributions p| . In particular, the discrete Wigner 
function has attracted a good deal of attention. There 
is no unique way of defining such a function: one family 
of methods maps the states of the Hilbert space onto a 
phase-space grid of 2N x 2N points (implying a redun- 
dancy of information) 0, S 13 j whereas it is also possible 
to map the states onto a (nonredundant) N x N phase 
space 0, 0, 0, 0, 0, Q . In this paper we shall inves- 
tigate the latter type of Wigner functions. 

Even-dimensional Hilbert spaces are of special inter- 
est to visualize the effects of quantum information pro- 
cessing. The ubiquitous qubit lives in a two-dimensional 
Hilbert space, and n qubits hence span a 2™-dimensional 
space. Wigner functions for two and three qubits have 
found applications in providing solutions t o, e. g ., state 
reconstruction |15|. quantum teleportation |l6lll7| . quan- 
tum optics [T3| . quantum computing |is|. and the mean 
king problem[2(J HJ |H |H H E| . 

In several of these contexts, mutually unbiased bases 
(MUBs) are also of interest |27| . We recall that each 
vector in one of these bases is an equal-magnitude su- 
perposition of all the vectors in any of the other bases. 
MUBs are central to quantum tomography and state re- 
construction 28], but are also valuable for quantum key 
distribution [IT |M EJ 111 . It is known that when the 



Hilbert space dimension N is a prime, or a power of a 
prime, there exist exactly N + 1 MUBs [33| . In this pa- 
per we shall mainly discuss the three-qubit case, for it 
is the smallest space consisting of qubits where differ- 
ent MUB structures, with respect to their entanglement 
properties, exist More precisely, in this space 

there are four different MUB structures corresponding to 
s = 0, 1, 2, and 3, respectively, where s is the number of 
triseparable MUBs. The fact that different structures are 
possible has strong implications for the physical imple- 
mentation of quantum tomographic measurements and 
quantum information protocols. It also affects the man- 
ner in which we can map states onto a phase space while 
retaining some properties that we consider "natural" in 
the continuous case. 



The rest of the paper is organized as follows: In Sec- 
tion [n] the continuous Wigner function and specifically 
the concept of translational covariance is recalled. In 
Section II I II a method for determining MUBs for three 
qubits is briefly reviewed. The concept involves operator 
sets based on primitive, one-qubit operations. In Sec. II VI 
we devise a construction of the Wigner function by asso- 
ciating the primitive operators (that generate a certain 
MUB) with translations in phase space. We show that 
this MUB has a unique structure, once we disregard the 
somewhat trivial degrees of freedom related to qubit la- 
beling and to association of lines or curves in phase space 
to states. In Section we derive two different Wigner 
functions from a different MUB structure, while in Sec- 
tion ^] we discuss the remaining two MUB structures, 
having s — or s — 1 triseparable bases. In this case, 
the principle of translational invariance no longer works 
in association with translations corresponding to single 
qubit MUB-generating operators. Finally, we make some 
concluding remarks in Section IVIII 
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II. THE CONTINUOUS WIGNER FUNCTIONS 

The continuous Wigner function has three properties 
that we are particularly keen on retaining when con- 
structing a discrete counterpart. The function is real, 
when integrated along any direction in phase space it 
yields a nonnegative function with unit area, i.e., a 
marginal probability distribution, and finally it is trans- 
lationally covariant. The latter property can be stated 
in mathematical terms as follows. Let W(q,p) be the 
Wigner function (expressed in terms of the position q 
and the momentum p) corresponding to a density matrix 
p, and let p' obtained from p by a displacement (qo,po) 
in phase space: 

p = exp[i(q p - p q)/h] p exp[-i(q p - poq)/h], (2.1) 

where q and p are the position and momentum operators, 
respectively. Then the Wigner function W correspond- 
ing to p' is obtained from W via the transformation 

W'(q,p)=W(q-q ,p-p ). (2.2) 

In other words, when the density matrix is translated, 
the Wigner function follows along rigidly. 

III. MUB STRUCTURES FOR THREE QUBITS 

Mutually unbiased bases can be constructed using a 
number of methods that depend on the dimensional- 
ity of the space. The main dividing lines are if the 
dimension is prime, a product of primes, or a power 

of a prime, and, in the latter case, if it is odd or 

even jMEaHlMliSlillillillillil- The prob- 
lem appears to be closely related to mutually orthogonal 
Latin squares 0, 23 an d to the existence of finite pro- 
jective planes of certain orders 0, 0] . As stated in the 
Introduction, we confine our study to the case of three 
qubits, that is, to an eight-dimensional Hilbert space. In 
this Hilbert space there exist four MUB structures, where 
the word "structure" denotes sets of MUBs where the ba- 
sis vectors are either triseparable (i.e., factorizable in a 
tensor product of three individual qubit states), separa- 
ble in one qubit and one maximally entan gled two-qubit 
state (i.e., biseparable) , or nonseparable |34 l35l| . The 
four structures are (3,0,6), (2,3,4), (1,6,2), and (0,9,0), 
where the labels indicate the number of triseparable, 
biseparable, and nonseparable bases, respectively. Note 
that all of them are of the form (s, b, 9 — s — b), where b 
is the number of biseparable bases. We shall primarily 
concentrate on the two first structures, since they corre- 
spond to MUBs where at least two bases are fully separa- 
ble so we can associate theses bases with local properties 
of the three qubits. This simplifies the analogy with the 
continuous Wigner function. 

We follow the construction algorithm worked out by 
Klimov, Sanchez-Soto, and de Guise in Ref. 22 to gener- 
ate, in a systematic way, a set of operators whose eigen- 
vectors constitute a MUB. If we take a spin-1/2 system 



TABLE I: Nine sets of operators defining a (2,3,4) MUB. 
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as our model for a qubit (described, apart from a factor 
h/2, by the Pauli operators <r), for three qubits it suffices 

to consider tensor products of , &y \ where 

the superscript i = 1,2,3 labels the qubit, in such a way 
that the tensor product contains one operator acting on 

each qubit. For example, <S> is a legitimate 

operator in this respect, which we shall write from now 
as <j x la y . 

We wish now to construct MUBs where two of the 
bases are generated by individual qubit rotations around 
two orthogonal axes, which we arbitrarily take as the x 
and the z axes. Hence, we define the MUB using the 
operators {ct z 11, ta z t, tla z } and {era,!!, 1(7x1) IIctz}. 
It can be shown that fundamentally, there is only a single 
way of accomplishing this, namely from the construction 
in Table |U 

The algorithm guarantees that the simultaneous eigen- 
states of the operators in each of the nine rows give a 
complete basis, and each basis is mutually unbiased to 
each other. From earlier work [3113 > we know that with 
respect to separability, Table |l] defines a (2,3,4) structure. 

Moreover, the table is uniquely defined by the 2x3 
entries in the three first columns of the first two rows. 
Indeed this is the case, since the operators in the first 
two rows are related by O r _ c = O r .c-20?- iC _3, where 
r = 1, 2 and c = 1, . . . , 7 denote the row and column 
of the operator, respectively, and must be taken modulo 
seven. Subsequently, the other rows can be expressed 
as O r . c — 02,cOi yC + r -3 for 2 < r < 9 (where all phase 
factors arising from the products have been neglected). 
In what follows, we write explicitly only the first 2x3 
elements of each basis to save space. 

We observe that if the first and second row are inter- 
changed, the new table defines the same structure. Our 
choice of axes is of course arbitrary, we could just as well 
denote the z axis by x and vice versa. If we permute the 
labeling of the axes, the result is the first set in Table ITU 
This looks different from Table HI and for fixed x and z 
axes it generates a different MUB than Table [TJ However, 
at the fundamental level, both are the same structure. If 
a certain election of the axes orientation is made, the 
MUBs are identical. 

Another arbitrary choice we have made is the labeling 
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TABLE II: The five different sets resulting from a rearranging of the three qubits of Table Q as 132, 213, 312, 321, and 231, 
respectively. Only the relevant 2x3 operators are shown. 



<r 2 ii ia z i iia z 


ia z i Ha z a z H 


ia z t a z tt lla z 




a z H ia z l 


Ha z ia z i a z lt 


a x ii iia x i& x i 


i& x i & x ii iia x 


M x l iia x & x ii 


ita x 


ta x i a x ii 


ita x a x ii ia x i 



of the qubits. If the qubits are numbered 1, 2, and 3, 
from left to right in Tabled rearranging the qubits (and 
the corresponding operators) in the orders 132, 213, 312, 
321, and 231 will then result in the five sets shown in 
Table ITT1 We see that, due to the simple structure of the 
tables, a permutation of the axes can also be seen as a 
relabeling of the qubits. 

The sets in Tables [I] and [H] are the only three- 
qubit, MUB generating matrices containing in one 
row the single qubit rotations around one axis, e.g., 
{(TjlLl, icJjlL, il<T z } (in any order), and in the next row 
the single qubit rotations around an orthogonal axis, e.g., 
{<7 X 11, ia x i, tta x }. Hence, with this restriction, only 
one MUB structure is allowed (up to a relabeling of the 
qubits). This is our first important result. 



IV. A THREE-QUBIT WIGNER FUNCTION 

When the Hilbert space is finite, there are several pos- 
sibilities for defining a Wigner function, as discussed in 
the Introduction. We follow the approach of Wootters, 
Gibbons, and Hoffman [H [3 and try to incorporate 
the three features discussed for the continuous case. To 
generate Wigner functions from the MUB in Table [I] we 
resort to elementary notions of finite fields. The field 
Fg with exactly eight elements can be seen as the set 
{0, 1, /i, . . . , /j?} (which is also an additive group), where 
the primitive element /i is a root of the following irre- 
ducible polynomial on Z2 (the integers modulo 2) 

8 3 + 6 + 1 = 0. (4.1) 

With this arithmetic, we have 

1 = n 3 + n 5 +v 6 , fi = fi 5 + ^ 6 , 

(4.2) 

/i 2 =/i 3 + ^ 5 , // = //+//. 

Note that the subset {fi 3 , fi 5 , fi 6 } defines a self-dual 
basis: 

tr(/V)=<%, i,je 3,5,6. (4.3) 

Here the trace of an element 9 in this field is defined as 

tr6 = 6 + 6 2 + 9 4 . (4.4) 

Following Refs. and 0, we subsequently associate 
the self-dual basis elements, in a two-dimensional phase- 
space coordinate representation, with the two groups of 



three operators, each one representing an "axis" in phase 
space: 

(^ 3 ,0) «-> a x tt, 

(^ 5 ,0) «-» la x l, (4.5) 

(m 6 ,o) «-> iis- X) 

and 

(0,/i 3 ) «-» &M, 

(0,/i 5 ) *-> ia-,1, (4.6) 

(0,/i 6 ) «-» Ua z . 

We map the eight-dimensional Hilbert space onto a 
8x8 discrete phase space, and label the x (horizon- 
tal) and z axes (vertical) with the field element sequence 
{0, 1, fx, [i 2 , . . . , fi 6 }. This gives us a phase-space coordi- 
nate in terms of field-element powers: each phase-space 
point has a unique coordinate, such as (/i 3 ,/x 6 ), as can 
be seen in Fig. ^ 

One can now construct lines and striations (sets of 
eight parallel lines, i.e., with no point in common) based 
on the algorithm for the MUBs. The striations are 
uniquely defined by (any) one line and the requirement of 
translational covariance. We shall, for definiteness, take 
the striation-defining line as the one that includes the 
origin, that is, a ray. Subsequently, the seven other lines 
belonging to each striation can be generated by trans- 
lating the ray in a manner that will be described below. 
In order to find a correspondence between the phase- 
space lines, the translations, and MUBs, we first define 
the unitary qubit-flip operators corresponding to the two 
sets of generating operators (that is, the set of opera- 
tors in a row in the respective defining table). The uni- 
tary operator exp(i7TCT x /2) ®i®i = C r 2; ®i(8)l (with 
U x = ia x ) rotates the first spin qubit around the x axis 
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FIG. 1: The striation-generating curves corresponding to the 
MUB construction defined by Table 1. 
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by the angle it. That is, the operator is an eigenoperator 
of the set {a^ll, lo^l, llo^}, and flips the first qubit of 
any eigenstate of {<r z ii, i<r z i, llcr z } into the orthogonal 
qubit state. In the same manner, local unitary flip oper- 
ators corresponding to the remaining five basis-defining 
operators {la x t,lla x ,a z ll,la z l,tl.a z } can be con- 
structed. 

In the phase-space representation these flip operators 
are represented by a translation by the corresponding 
field element "vectors", e.g., U x <g> 1 <g> 1 *-> translation 
by addition of (/x 3 , 0) , and 1 ® U z <8> 1 <-> translation 
by addition of (0,/x 5 ). The ray associated with the set 
{o^ll, ia x i, iii7 x } is hence generated by translation of 
the origin {0, 0} by all additive (modulo 2) combinations 
of (^ 3 ,0), (n 5 ,0), and (^ 6 ,0). It is clear that this will 
result in the horizontal line (ray) of Fig. (The origin 
is marked with an o because it belongs to all striation- 
generating rays.) It is also evident that if every point 
on this ray is transformed by any combination of the 
three translations, the ray will remain invariant. This 
reflects the fact that the state associated with the ray 
is an eigenvector of the corresponding unitary qubit flip 
operators. 

The rest of the lines of this striation can be obtained 
by translation by all the combinations of the three vec- 
tors (0,/x 3 ), (0, fj, 5 ), and (0,/x 6 ) (respecting the modulo 2 
arithmetic and the field algebra). This results in a set 
of eight parallel lines - the striation corresponding to the 
basis generated by the operators {(Txll, Ict^I, IIct^}. In 
the same manner, one can construct the striation (con- 
sisting of straight vertical lines) representing the MUB 
generated by {<r z ll, 1<t z 1, 11<t z }. The ray belonging to 
this set is labeled with 2s in Fig. ^ 

To generate the third striation-generating "ray" , we 
follow the recipe given above [42j and simply multiply the 
first two rows of Table [I] columnwise. In doing so, we ig- 
nore overall phase factors. We get the third set of MUB- 
generating operators {a y tt,ta x a z ,ta z a x }. These op- 
erators correspond to the phase-space operators (/x 3 , fi 3 ), 
(/x 5 , /j 6 ), and (/i 6 ,^ 5 ). The "ray" generated by displace- 
ment of the origin by all the combinations of these trans- 
lations is labeled with 3s in Fig. 2] 

The reason we have used the word "ray" within quo- 
tation marks in the paragraph above is that "ray" 3 is in 
fact a curve in phase-space. If a point belonging to it is 
denoted (a, then the curve is parametrically defined 
as 

a = [i 3 k + [i 5 k 2 + k 4 , (3 = fi 2 n + n 2 + ^ 4 k 4 , (4.7) 

where k is a parameter running over the field elements. 
This parametrization corresponds to the explicit expres- 
sion 

f3 2 + fi/3 = a 2 + fjux, (4.8) 

or, equivalently, 

(3 = fjpa + //a 4 + (j, 3 a 2 , (4.9) 



from which it is clear that a and (3 are not linearly re- 
lated. In the following we will refer to a curve passing 
through the origin as a homogeneous curve. 

From the table we see that, as required by a set of 
MUBs, each curve crosses the curves of the other stria- 
tions only once, at the origin in Fig. ^ by our choice of 
depicting only homogeneous curves. Every full striation 
can be generated from the homogeneous curve by trans- 
lation in the horizontal direction by adding, in sequence, 
all the combinations of (/i 3 ,0), (/i 5 ,0), and (/x 6 ,0). In 
physical space, this corresponds to all possible combi- 
nations of spin flips around the x axes. The striations 
can equally well be generated by vertical translation by 
all the combinations of (0,^i 3 ), (0, fj, 5 ), and (0,/x 6 ), cor- 
responding to spin flips around the z axes. In fact, we 
can generate any striation from the corresponding homo- 
geneous curve by the translations corresponding to any 
of the other curve-generating operators. For example, 
the operators (/x 3 ,/Lt 3 + /i 5 ), (/x 5 ,// 3 ), and (/x 6 ,/x 6 ), corre- 
sponding to the first three columns of the bottom row of 
Table U will also generate the same striation. This follows 
from the fact that the relation between the bases is the 
same for any two of them. The fully separable bases are 
not particular in this respect. 

To obtain a Wigner function from the striations one 
more step is needed, namely, to associate each line in 
phase space with a state. Again this involves an arbi- 
trary choice that does not lead to anything fundamen- 
tally new. As described in Ref. ITil we can associate 
any basis eigenstate to any line in a striation, but once 
this choice is made, the state associated to a line ob- 
tained by a particular translation must correspond to 
the state obtained from the first by the corresponding 
spin-flip transformation. For example, if we associate 
the vertical ray in Fig. (labeled with 2s) with the 
state | TTt)z ( m a spin- 1/2 representation with z as our 
spin axis), then the vertical line that includes the point 
(ji 3 + [i 6 , 0) = (^ 4 , 0) U x <g> 1 <8> U x represents the state 
U x <8> i ® U x \ TTT) 2 = I IU)*< As there are N + 1 stri- 
ations, and for each striation we can associate the gen- 
erating ray with any of iV state vectors, there are N N+1 
quantum nets associated with this choice of phase-space 
bases (single qubit rotations around the x and z axes). Of 
these, one can group the possibilities in N N_1 equivalence 
classes each containing N 2 elements 14]. The elements 
are simply the N xN possible discrete translations of the 
space using all combinations of the horizontal and verti- 
cal translation operators. Nevertheless, we stress that in 
our approach, in contrast to that of Wootters [l3l[T4| . we 
associate curves in phase space (and not merely lines) to 
states. 

With the prerequisites above, for any phase-space 
point (a, (3) we can define Hermitian phase-space point- 
operators A( Q>| g) as 

9 

i(a,/3)=^/3fc,(a,/3)-i (4-10) 

fc=l 
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where pk,(a,0) 1S the density matrix associated with the 
curve belonging to striation k passing through (a,/3). 
From the mutual unbiasedness of the orthonormal bases, 
it follows that 

Tr[i (Qi/3) ] = 1, 

(4.11) 

where Tr (with upper case) denotes the ordinary Hilbert- 
space trace operation. The discrete Wigner function 
W(a, (3) of the density matrix p is then defined as 

W(a,(3) = ±Tr[pA (a , /3) ]. (4.12) 

This leads to the relations 

5^W(o,/3) = l, 

a, 13 

p = J2w(a,/3)A (a>/}) , (4.13) 

a, [3 

Tr(pp')=8^2W(a,p)W'(a,0), 

a, 

where W'(a, (3) is the Wigner function of the (arbitrary) 
density matrix p'. It is seen that the phase-space point- 
operators, and consequently the Wigner function, fol- 
low naturally from the construction algorithm outlined 
above, based only on the generating operator tables, the 
irreducible polynomial (|4.1|l . and the (arbitrary) associ- 
ation between basis states and phase-space rays. 

The phase-space structure corresponding to the first 
set in Table ITT1 is obtained by making a mirror image of 
Fig. with respect to the diagonal through the origin. 
As noted above, this stems from the fact that the re- 
placement x <-> z generates this set in Table ITT1 from the 
one m Tabled and vice versa. The phase-space struc- 
ture corresponding to the other sets in Table |H] can be 
derived in the same manner. However, the interrelation 
between the structures through qubit permutations is not 
as easy to unveil in phase space as in physical space, as 
demonstrated by Fig. [21 



V. THE (3,0,6) MUB-STRUCTURE WIGNER 
FUNCTION 

As we have seen above, individual rotations of the 
qubits around the x and z axes only generate one MUB 
structure. To proceed further, one needs to consider si- 
multaneous rotations of two or more qubits. The simplest 
one, in which we allow simultaneous rotation of the first 
and the second qubits around the x and z axes, is shown 
in Table ITTT1 

TABLE III: A table defining a (3,0,6) MUB structure. 

0" Z O" Z 1 11<7 Z 0" Z 1]1 



The rows of the table have the same eigenvectors as 
the corresponding rows of Tables □ and QH above. Mul- 
tiplying the rows columnwise, we obtain the operators 
{(TytTyl, liiTy, CTj,!!}. These operators have the associ- 
ated phase-space translations (p 3 +p 5 , p 3 + p 5 ), (p 6 ,p 6 ), 
and (p 3 ,p 3 ). It is evident that all the three operators 
make "diagonal" translations in phase space. It is also 
evident that the corresponding MUB is triseparable: it 
is a (3,0,6) structure. The striation-generating rays of 
the phase-state structure are displayed in Fig. |3J 

This table has a particularly simple structure although 
it is built using the same algorithm as all the preced- 
ing phase-space structures. The rays in Fig. |3 left side, 
are really rays, and hence all striations will consist of 
lines. That is, the different rays are generated by the 
equations f3 = Xa and a = 0, where the "slope" A for 
each ray is a field element. Moreover, the rays are "visu- 
ally straight". This, e.g., means that Table ITTT1 generates 
the same MUB under the interchange of axes x <-> z. 
We can once more generate six "different" tables, rep- 
resenting the same physical structure, by relabeling the 
qubits. The simplicity of this figure is a chimera, how- 
ever, because only the rays are diagonal and "visually 
straight". Translating the rays by, e.g., flipping the first 
qubit around the z axis (corresponding to a translation 
by (0, p 3 ) yields the lines in the right part of Fig. [3] In 
this set of lines, one from each striation, only the hori- 



2 


6 


5 


8 


3 


7 


4 


9 


2 


5 


4 


7 


9 


6 


3 


8 


2 


9 


8 


4 


6 


3 


7 


5 


2 


7 


6 


9 


4 


8 


5 


3 


2 


4 


3 


6 


8 


5 


9 


7 


2 


8 


7 


3 


5 


9 


6 


4 


2 


3 


9 


5 


7 


4 


8 


6 


o 


1 


1 


1 


1 


1 


1 


1 



2 


5 


7 


6 


8 


4 


9 


3 


2 


7 


9 


8 


3 


6 


4 


5 


2 


8 


3 


9 


4 


7 


5 


6 


2 


6 


8 


7 


9 


3 


1 


2 


2 


9 


4 


3 


5 


8 


6 


7 


2 


4 


6 


5 


7 


3 


8 


9 


2 


3 


5 


4 


6 


9 


7 


8 


o 


1 


1 


1 


1 


1 


1 


1 



FIG. 2: The striation-generating lines corresponding to the 
MUB construction defined by the second and fourth sets in 
Table 2. 
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FIG. 3: The striation-generating rays corresponding to the 
MUB construction defined by Table 3, left, and the lines re- 
sulting from displacement of the rays by (0, p 3 ), right. 
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TABLE IV: Another table denning a (3,0,6) MUB structure. 
a z W a z a z \ Wa z 

G X W & x a x 1L UlTa, 
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FIG. 4: Homogeneous curves corresponding to the MUB con- 
struction defined by Table 4. 

zontal and the vertical lines remain "visually straight" . 

The interesting point is that there exist a physically 
different (3, 0, 6) MUB structure, also based on two single 
qubit-rotations and one two-qubit rotation. It is shown 
in Table ITVI 

The homogeneous curves corresponding to this struc- 
ture are displayed in Fig. A significant result is that 
the structures represented by Tables ITTll and llVl are the 
only physically inequivalcnt ones that one can obtain 
from one double rotation and two single qubit rotations. 

VI. MUB STRUCTURES (1,6,2) AND (0,9,0) 

It is known that two more MUB structures exist in 
the three-qubit space, the (1,6,2) and (0,9,0) [M H||. 
The first has only one triseparable basis, whereas the 
second has none. It is still possible to associate these with 
Wigner functions in a similar manner, the only difference 
is that the corresponding Wigner functions will have only 
one or no "axes" (striations composed of lines). This 
means that most striations will consist of curves with 
two, and sometimes four, points on the same horizontal or 
vertical coordinate. To demonstrate this we use a (1, 6, 2) 
structure, which can be generated from the operator sets 
in Table El 

We see that the set involves "simultaneous" rotations 
of the first and middle qubits around the x and z axes, 
effectively resulting in a rotation around the y axis. In 
phase space this means that by choosing single qubit ro- 
tations around the x and z axes as our "primitive" oper- 
ations, no set of striation-generating translations decou- 
ple in horizontal and vertical translations. Therefore, the 
correspondence with the translation and boost operators 
q and p, respectively, in the continuous case is lost. The 
homogeneous curves corresponding to this structure are 
depicted in Fig. [3] We see that since there is no trisep- 
arable basis involving only rotations around the x axis, 
there is no horizontal ray. Instead, the curve labeled 1 



TABLE V: A table defining a (1,6,2) MUB structure. 

1<t z 1 a z a z \ Wa z 
b x o y \ Wg x a y b z \ 
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FIG. 5: Homogeneous curves corresponding to the MUB con- 
struction defined by Table 5. 

can parametrically be defined as 

a = /i 2 K 4 , (3 = fi 2 n + k 2 + /in 4 , (6-1) 
which corresponds to the explicit expression 

P 2 + ^(3 = fi 3 a 2 + fi 2 a, (6.2) 
or, in an equivalent form, 

13 = + /j 5 a 2 + //a 4 . (6.3) 

This in turn implies that the other curves have several 
points on the same horizontal line. For example, curve 6 
has four points on the horizontal lines with z-coordinatcs 
and /i 4 . This reflects the fact that this is a curve invari- 
ant under the two horizontal translations (/i 3 + /x 6 , 0) and 
(/i 5 ,0). By applying the translation operators (0,/i 5 ), 
(0, /i 3 + [i 5 ), (0, /x 6 ) one can subsequently generate stria- 
tions from curves 1 and 3-9. However, this set of trans- 
lation operators leave ray 2 invariant. To generate the 
whole striation from ray 2 we can, e.g., use the transla- 
tion set (fi 3 + (U 5 ,// 5 ), (aAO), and (/i 3 ,^ 3 + A* 5 ), which 
are the generating translations of curve 1. 

To finally depict a "severely disordered" three-qubit 
Wigner function (at least to the human eye) , we use Table 
IVII generating a (0,9,0) MUB, which gives rise to the 
nine striation generating homogeneous curves depicted in 
Fig. In spite of its disordered appearance, the Wigner 
function defined from Table IVII inherits all the desirable 
features from the continuous Wigner function except for 
(visually) straight axes corresponding to q and p. This 
follows from the fact that all the MUB have entangled 
basis vectors. 

TABLE VI: A table defining a (0,9,0) MUB structure. 

(Tajll \tj z tj y TLa y a z 
(Tyll ta x a z ta z a x 
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FIG. 6: Homogeneous lines corresponding to the MUB con- 
struction denned by Table 6. 



VII. CONCLUDING REMARKS 

We have discussed the relation between MUBs hav- 
ing different structures with the respect of the entan- 
glement properties and translationally covariant Wigner 
functions for three qubits. We have shown that there 
exist three fundamentally different constructions of the 
discrete Wigner functions if single qubit rotations around 
two orthogonal axes are used as the generating opera- 
tions. To construct the Wigner function we need only 
a MUB generating table and a field-element generating 
irreducible polynomial. Other constructions can then 
be generated by qubit permutation and by associating 
curves in the striations with different states (the choice 
being arbitrary). The number three is surprisingly small, 
and shows that requiring translational covariance im- 
poses severe restrictions. 

We have also shown that for three-qubit MUBs that 
have only one or none triseparable bases, the method 
based on the above-mentioned generating operations still 
works. Here we do not know how many fundamentally 
different constructions exist. However, these MUBs are 
perfectly legitimate from a physical point of view, al- 
though they have a more complex appearance than the 
Wigner functions based on (2,3,4) and (3,0,6) MUBs. 

We have not addressed here the interesting question 
of the factorizability of the Wigner function, that has 
been previously considered by Durt [5(J and Pittenger 
and Rubin |5l|. This is certainly relevant in the frame- 
work of quantum tomography, particularly from an ex- 
perimental viewpoint. In principle, it would suffice that 
one phase-space point is trifactorizable in order that the 
other operators are, since the translations themselves are 
factorizable in virtue of translational covariance. Unfor- 



tunately, all our efforts to determine if there exist some 
triseparable phase-space point operator have failed. To 
the best of our knowledge, the problem thus still remains 
open. 

In general, the curves that define a MUB on F 2 3 form 
one-dimensional Abelian structures, which can be conve- 
niently parametrized as 

a(n) = V1K + V2K 2 + v?,k A , (3(n) = riiK + r]2K 2 + rj 3 n 4 , 

where k is a parameter and the coefficients Vj and rjj take 
values on F 2 3 , so that 



o{k + k) = ci(k) + q;(k), 



/8(k + k') = P{k)+/3(k'). 



(7.2) 



Hence, such curves are the simplest generalization of 
"straight" Abelian structures (rays) of the form 

a = 0, or = Xa, (7.3) 

or equivalcntly 

ol{k) — r)K, P(k) = Qk, (7-4) 

where r\ and £ are fixed field elements. 

We finally observe that the structures studied byWoot- 
ters llalbl (and also by Bandhyopadhyay et al [4l[ and 
Durt |I(|) always assume two trifactorizable bases, so 
that only the cases (2,3,4) and (3,0,6) are possible. The 
existence of a third triseparable basis depends on the 
choice for the field basis: in the self-dual used in this pa- 
per, this is always the case [so we are lead automatically 
to the (3,0,6) structure], while other choices bring the 
(2,3,4). In this respect, it is interesting to note that the 
(3,0,6) is the only three-qubit MUB phase-space struc- 
ture (depicted in Fig.|3J) consisting solely of straight lines. 
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